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Introduction 

Pyrometer is a favorite method to do remote temperature measurement in research and development. 
One-color, two-color and the disappearing filament pyrometers*" are most common, multicolor and 
multiwavelength pyrometers are being introduced recently* 2 "". All these pyrometers invariably require, 
in one form or another, information concerning emissivity, the medium transmissivity, their ratio at 
some two spectral regions, the instrument’s calibration constant, etc. for their operation. This 
information can come from hand books, from the manufacturer or in some instances, from results of 
dedicated separate experiments. Often this information is sample or instrument specific. Sometimes 
this information, though obtained from a special experiment, is obtained using a separate sample 
rather than the one pursued for temperature measurement. Then, there would be the question of 
variability from sample to sample and variability from batch to batch in the sample used. Also, 
previously determined calibrations can change with time, and the only way to reduce uncertainty is 
to perform the calibration more frequently or even immediately before the experiment. We have 
developed a multiwavelength pyrometer, which eliminates the need to supply the necessary emissivity 
and/or transmissivity information and the instrument calibration constants ahead of time. The 
pyrometer calibrates itself from its first cycle data. 


Theory 

Planck’s law of black body radiation (Eqn 1) is the foundation of pyrometry, with c 1; a, the radiation 
constants, X the wavelength and T the black body temperature. 
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In its most basic form, a pyrometer is a spectrometer which records a voltage V(X.,t) at time t in 
response to radiation emitting from a black body. In general, the radiation source is not a black body, 
and the output voltage V(k,t) depends on the spectral emissivity of the emitting body, x x the 
transmissivity of the optical medium between the emitting body and the detector, and the instrument 
constant g x given by Eqn 2 
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The temperature is assumed to be a function of t. Ideally, if g, £ and x are known at wavelength X, 
temperature measurement is reduced to measurement of V(k,T) only. In general, g, e and x are not 
known We show that though temperatures TXt^Tit^ .. *T(t,)* ~*T(t n ) are unknown, spectra of 
wavelengths X v X 2 , ..X,.. V N>2, obtained at times t v t*, .. t,, .. t n , contain enough information to 
determine everything that is required in pyrometry. Consider at time t, 2 wavelengths, label them 
and their associated quantities by Xr and ^ 
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Solve for T(t), replace subscripts X R , X i by just R and i in g, £ and x, we have 
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Elimin ating temperature T(t) yields Eqn 6 from Eqn 5. Letting AR=gR£ R x R and Aj-gjEjXj, Eqn 6 becomes 
Eqn 7, this is the equation which we will use. 
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Assume A R =gR£ R x R is constant, i.e. time and temperature independent, the A i =g i £ i x i s may dependent 
on temperature, then a plot of the quantity on the left hand side vs the quantity inside the curly 
brackets on the right hand side of Eqn 7 will produce a straight line of slope l/A i =I/g i e i x i and intercept 
. l/A i =- 1/gjEjXj . Do this for all &R. The spectral quantity Aj is a function of Ar. Referring to Eqn 5 or 
its equivalent, Eqn 8, use the generated Aj to calculate the temperature T(A.,t) at time t. 
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In theory, these calculated temperatures should all equal to each other and be independent of 
wavelength. In actual data, there are variations, an average T(t) is obtained according to Eqn. 9. 
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This is done for all the spectra (i.e all t) to obtain average temperatures as determined by a particular 
value of Ar. The correct value A R is determined by a least squares procedure as follows: 


1) 
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7) 


Choose a value for Ar, 

Plot the data according to Eqn 7 to determine the Aj from the slopes. 

Use the so determined Aj to calculate the temperatures T(t) according to Eqns. 8 & 9. 
Transform the spectra V(X,,t) into a large data set (x,y). The wavelength \ is transformed into 
c/XTit), the transformed wavelength, the voltage V(^,t) is divided by Aj, and then by T(t) 5 , the 
5th power of the spectrum temperature, according to the prescription 
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The transformed (x,y) data obey the generalized non-dimensional Planck function. The (x,y) 
data is fitted to the Planck function of Eqn. 10 by calculating the residual. The residual 2, is 
defined as the sum of the squares of the difference between the transformed y, and the 
calculated y evaluated by substituting the transformed Xj in the y equation in Eqn 10 for all 
the data. 

A new value for Ar is selected and steps (2 to 5) repeated. 

The particular value of A R that produced the least £ is the correct one. 


Thus in one step, A R , also A; and T(t) are determined. In one stroke, everything that is needed in 
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pyrometry for temperature measurement is determined. Temperatures at any time in the past or in 
the future are determined from each A, according to Eqn 4 or its equivalent Eqn 8. As is evident, the 
pyrometer requires no prior calibration. When N=2, A; is defined in terms of A R and the slope of the 
experimental data. The correct value of A R is similarly determined by least squares curve fitting the 
N=2 smaller transformed (x,y) data sets at different temperatures exactly as described above to Eqn 
10. Greater redundancy is achieved when N>2. 

Experiment and Results 

Two sets of results are presented. The first is temperature measurement of a black body furnace 
viewed through a window of unknown transmisivity, the spectrometer has not been calibrated. The 
second one is the temperature measurement of a zirconia ceramics of unknown emissivity, the 
spectrometer has been calibrated. 

(I) Unknown Transmissivity Case 

The experimental arrangement is shown in fig. 1. This is the case of unknown g x , unknown t x and 
^=1. A commercial black body furnace operated between 1000 °C and 2000 °C is used. Its graphite 
heating element and cavity are protected against oxidation by an inert argon gas curtain. When 
viewing the black body cavity directly is not required, a calcium fluoride window is normally positioned 
at the furnace opening to reduce the argon gas consumption. In the described experiment, the calcium 
fluoride window was always left in place. The black body furnace temperature was increased in steps 
by gradually changing the controller from its lowest to its maximum temperature set point. An 
indicator is present to show whether the pre-set temperature has been reached. When that occurred, 
a spectrum was recorded. The spectrum spanned the spectral region from 0.6 to 4.5 pm. It recorded 
the direct voltage output of the indium antimonide detector. 

Results Six of these spectra are shown in fig. 2. The spectrum minim a are due to atmospheric 

C0 2 and HjO absorptions in the optical path between the detector and the black body source. 
Following the analysis above, Aj, is chosen to be at 2 pm. An arbitrary initial value for A R is used, 
plots according to Eqn 7 are made. Fig. 3 shows the case for V=2 pm, and A, =1.5 pm, A R =0.1534. It 
is indeed a straight line. The slope and intercepts of plots like this at other wavelengths are obtamed 
using least squares method. They are plotted in fig. 4. Individual Aj is obtained from the slope data 
because less noise is present, the result is shown in fig. 5. The Ai are now used to calculate the 
temperature of each spectrum according to Eqn 8 and shown in fig. 6. These calculated temperatures 
are almost independent of wavelength. The explosion at the shortest wavelength is due to poor signal 
to noise there, and we excluded data at wavelength shorted than 0.8 pm in the analysis. An average 
is obtained for each spectrum, i.e. at time t. When the value of changes, the slopes also change, 
hence A,, which depend on the slopes, also change and finally the temperatures at each wavelength, 
calculated using Aj and the spectral data, also change. Of the many possible values that A R can 
assume, we identify the correct one using the least squares curve fitting method, when the reduced 
(x,y) data set is used to fit the non-dimensional Planck function of Eqn 10. The least squares condition 
is obtained when A K =0.1534 and the result is shown in fig. 7. Six spectra were used. More spectra 
could be used. Using more spectra would increased the (x,y) data set size, and greater redundancy, 
leading to better accuracy. The maximum number of spectra that we can use is limited only by 
computer memory. With the determination of just one parameter A R , the Aj at all other wavelengths 
as weU as the temperatures of all the spectra that were included in the analysis are also determined 
at once. Once A, is determined, the temperature at any time in the past or in the future are 
determined immediately from the measured spectra. Alternately, this procedure can be repeated as 
often as required to update the calibration. 

The measured temperatures are 2234, 2143, 1965, 1804, 1667, 1550, 1460, 1379, 1313 and 1256 K, 
when the black body furnace controller pre-set temperatures are 2293, 2133, 1955, 1793, 1668, 1574, 
1448, 1378, 1318 and 1268 K. The results are shown in fig. 8. The agreement is within 0.5 %, except 
at the pre-set temperature of 2293 K, when the measured temperature is 2.5 % lower. We believe that 
the measured temperature is actually more correct, because during the experiment, according to the 
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controller indicator, the black body furnace never reached its set point temperature. The controller 
required a power line voltage of 220 V, but the actual line voltage is only 208 V. The temperature 
2234 K was the eq uilib rium temperature when the controller was delivering its ma x i m u m power 
output under the experimental line voltage condition. 

For the N=2 case, we choose Xr= 2 pm, ^=1.5 pm. With 6 spectra, there are now 12 data points in the 
transformed data set (x,y). Least squares curve fitting these data to Eqn 10 gives the value 0.1471 
for A r . The fitted curve is shown in fig. 9, and the resulting temperatures for the spectra are 2262, 
2169, 1987, 1822, 1683, 1564, 1472, 1390, 1323 and 1265 K, which are plotted in fig. 10. The largest 
error in these temperatures is less than 1.7 % of the pre-set temperatures. The use of only two 
wavelengths greatly reduced the data acquisition time by a factor of 10 or even 100. 

(II) Unknown Emissivity Case 

This is the case of unknown e^, known g x and t x = 1. Emission spectra of a 250 pm thick zirconia (Zr0 2 ) 
thermal barrier coating (TBC) are measured using the pyrometer arrangement of fig. 11 from 1.3 to 
14.5 pm. They are shown in fig. 12. Zr0 2 TBCs have wavelength dependent emissivity <5) , being very 
email (£<0.4 or 0.3) at short wavelengths, increasing to almost unity at the longer (X>10 pm) 
wavelengths. Zr0 2 TBC temperature measurement using the traditional 1- and 2-color pyrometers 
are known to be difficult. The TBC is flame sprayed on a 25 mm diameter, 6 mm thick metal disk. 
A 12 mm deep small hole is drilled along a radius 3 mm from the disk surface into which a 250 pm 
type G thermocouple (TC) is positioned to provide reference temperatures. The disk is placed about 
25 mm inside the opening of a black body furnace opening whose temperature is carefully controlled. 
These spectra (fig. 12) were recorded when the imbedded TC indicated steady readings (at 1110, 1001, 
892, 791, 693 and 602 K). 

Results In this experiment, the pyrometer is calibrated in the sense that g x is known and the 
measured spectrum is in energy units. The Zr0 2 TBC data of fig. 12 are analyzed according to Eqn 5 
assuming t k =T;=1, the emissivity e, is still not known, they will be determined together with the 
unknown temperature of each spectrum. The 250 pm thick Zr0 2 TBC transmits only less than 1% 
radiation in a very restricted spectral region (fig. 13), the spectra in fig. 12 must be due to emission 
or reflection. Reflection is not important in this experiment. Xr is chosen to be 5 pm and the analysis 
is done exactly as in case (I) for the unknown transmissivity case. Six spectra were chosen, slopes 
were determined, an initial value for e R is chosen, the individual £j are calculated in terms of e R and 
the slopes. The ^ are used to evaluate an average temperature for each spectrum, and similar to case 
(I), the spectra data are transformed by dividing each spectral datum by £*, then by T(t) s , the 5th 
power of the temperature, also the wavelengths are transformed into Cj/XT(t), and the transformed 
data set (x,y) is least squares curve fitted to the invariant curve. The resultant slope and intercept, 
emissivity and fitted curve are shown in figs. 14, 15 and 16. The emissivity shown in fig. 15 is very 
cimilar to that reported by Liebert <6) . The temperatures were determined to be 1096, 999, 895, 793, 
693 and 605 K differing from the TC values by less than 2% (fig. 17). 

Conclusion 


The multiwavelength pyrometer successfully measured the temperatures of a black body furnace 
viewed through a transparent window of unknown transmissivity, and also measured the temperature 
of a zirconia ceramic of unknown emissivity. The measurement error is less than 0.5 % for the first 
case. This is significant because the pyrometer is not previously calibrated. In the latter case, the 
error is less than 2 %, this is significant because the measurement is done without prior knowledge 
of this non-gray (wavelength dependent emissivity) material in the region where the emissivity is 
inherently small and interference from other sources is always present. By choosing just 2 
wavelengths in a good spectral region, temperature determination can be done and is much simplified 
with a small increase in error. 
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Fig. 1 

Arrangement of Experiment. 


PYROMETER 





W I NDOW 



Fig. 2 

Voltage spectra of black body furnace viewed through calcium fluoride window. 
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Fig. 5 

Calibration constant A; as a function of wavelength. 



Fig. 6 

Temperatures calculated from spectral data and A; as a function of wavelength. 
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Fig. 7 

Least squares curve fitting of the transformed data set to the non-dimensional 
Planck function. 


Transformed Spectra Fitted 
to Generalized Planck Function 



Fig. 8 

Plot of pyrometer measured temperature vs black body pre-set temperature. 



9 



Fig. 9 

Least squares curve fitting of the transformed N=2 case data set to the non- 
dimensional Planck function. 
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Fig. 10 

Plot of pyrometry measured temperature vs black body preset temperature for the 
N=2 case. 
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SLOPE AND INTERCEPT 


Fig. 13 

Tr ansmis sion of different thickness zirconia TBCs. 



Fig. 14 

Slope and intercept results at different wavelengths. 
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EMISSIVITY 


Fig. 15 

Emissivity of zirconia TBC as a function of A,, e R =0.5. 
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Fig. 16 

Zr0 2 TBC spectra transformed into generalized Planck equation. 
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Fig. 17 

Pyrometry and thermocouple measured temperatures. 
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